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The Fokker Planck (FP) equation has been solved to study the interaction of non-equilibrated 
heavy quarks with the Quark Gluon Plasma (QGP) expected to be formed in heavy ion collisions at 
RHIC energies. The solutions of the FP equation have been convoluted with the relevant fragmen- 
tation functions to obtain the D and B meson spectra. The results are compared with experimental 
data measured by STAR collaboration. It is found that the present experimental data can not 
distinguish between the pr spectra obtained from the equilibrium and non-equilibrium charm dis- 
tributions. Data at lower pr may play a crucial role in making the distinction between the two. The 
nuclear suppression factor, Raa for non-photonic single electron spectra resulting from the semilep- 
tonic decays of hadrons containing heavy flavours have been evaluated using the present formalism. 
It is observed that the experimental data on nuclear suppression factor of the non-photonic electrons 
can be reproduced within this formalism by enhancing the pQCD cross sections by a factor of 2 
provided the expansion of the bulk matter is governed by the velocity of sound, c s ~ l/\/4. Ideal 
gas equation of state fails to reproduce the data even with the enhancement of the pQCD cross 
sections by a factor of 2. 

PACS numbers: 12.38.Mh,25.75.-q,24.85.+p,25.75.Nq 
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I. INTRODUCTION 



The nuclear collisions at Relativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider 
(LHC) energies are aimed at creating a phase of mat- 
ter where the properties of the matter is governed by 
quarks and gluons Such a phase of matter, com- 
posed of mainly light quarks and gluons - is called 
quark gluon plasma (QGP). The study of the bulk 
properties of QGP is a field of high contemporary 
interest and the heavy flavours namely, charm and 
bottom quarks play a crucial role in such studies, 
because they are produced in the early stage of the 
collisions, they are not part of the bulk properties 
of the system and their thermalization time scale is 
larger than the light quarks and gluons and hence 
can retain the interaction history more effectively. 

The successes of the relativistic hydrodynamical 
model 0, H| in describing the host of experimental 
results from RHIC [H indicate that the thermaliza- 
tion might have taken place in the system of quarks 
and gluons formed after the nuclear collisions. The 
strong final state interaction of high energy par- 
tons with the QGP i.e. the observed jet quench- 
ing 0, @ and the large elliptic flow (t> 2 ) @, § in 
Au+Au collisions at RHIC indicate the possibility 
of fast equilibration. On the one hand the experi- 
mental data indicate early thermalization time ^0.6 
fm/c Q on the other hand the pQCD based calcula- 
tions give a thermalization time ~ 2.5 fm/c (Toj(see 
also [ll|). The gap between these two time scales 
suggests that the non-perturbative effects play a cru- 
cial role in achieving thermalization. It has also been 
pointed out that the instabilities [HI, 0, [TJ, 
may derive the system towards faster equilibrium. 



The two pertinent issues regarding the equilibra- 
tion which will be addresses here are (i) do the the 
heavy quarks achieve equilibrium and (ii) in case 
they achieve equilibration, can the equilibrium be 
maintained during expansion of the system. The 
second issue will be addressed first. 

We make a rather strong assumptions that the 
heavy quarks produced initially is in thermal equi- 
librium and check whether it can maintain the equi- 
librium during the entire evolution processes by com- 
paring their scattering rates with the expansion rate 
of the matter. This issue will be addressed with 
different equation of states (EoS) which affects the 
expansion rate. In case the heavy quarks are un- 
able to maintain the equilibrium then the anal- 
ysis of the transverse momentum of the mesons 
carrying heavy flavours can not be done by us- 
ing thermal phase space distribution. The analy- 
sis will require non-equilibrium statistical mechani- 
cal treatment. We solve Fokker-Planck (FP) equa- 
tion [H [13, H E! [H M, EH, HI to address this 
issue as discussed below. 

The perturbative QCD (pQCD) calculations indi- 
cate that the heavy quark thermalization time, 
is larger than the light quarks and gluons ther- 
malization scale Ti. Gluons may thermalize before 
up and down quarks [HI, [24[ , in the present work we 
assume that the QGP is formed at time Tj. There- 
fore, the interaction of the non-equilibrated heavy 
quarks (Q) with the equilibrated QGP for the time 
interval Tj < r < can be treated within the am- 
bit of the FP equation i.e. the heavy quark can be 
thought of executing Brownian motion in the heat 
bath of QGP during the said interval of time. The 
solution of the FP equation can be used to study pt 
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spectra of heavy mesons in the spirit of blast wave 
method. 

In the next section we address the issues of ther- 
malization in a rapidly expanding system. The re- 
sults indicate that the heavy quark can not main- 
tain the equilibrium at RHIC and LHC energies dur- 
ing the entire evolution history of the QGP. This 
demands the treatment of the problem within the 
framework of non-equilibrium statistical mechanics, 
which is discussed in section III. Section IV is de- 
voted to summary and conclusions. 



II. THERM ALIZATION IN AN 
EXPANDING SYSTEM 



We consider a thermally equilibrated partonic sys- 
tem of quarks, anti-quarks and gluons produced 
in relativistic heavy ion collisions. We would like 
to study whether the system can maintain ther- 
mal equilibrium when it evolves in space and time. 
Relativistic hydrodynamics (with boost invariance 
along longitudinal direction and cylindrical symme- 
try) have been used to describe the space-time evo- 
lution. For this purpose, the scattering time scale 
(T sca tt) of the partons are compared with the expan- 
sion time scale {r cxp ). For maintenance of thermal 
equilibrium the following criteria should be satisfied: 



(1) 



where a ~ 0(1) is a constant. The criteria given in 
Eq. [3 is reverse to the one used to study the freeze- 
out of various species of particles during the evolu- 
tion of the early universe [25| (similar condition is 
used in (2(| in heavy ion collisions also). 

The T sca tt is determined for each partons by the 
expression 
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2J aijVijUj 



(2) 



where cry is the total cross section for particles i and 
j, Vij is the relative velocity between the particles i 
and j and rij is the density of the particle type j. 

To calculate the scattering time we use the fol- 
lowing processes gg -> gg, gg -> qq, q(q)g q(q)g, 
qq — > qq, qq — > qq for light flavours and gluons [271 ] . 
Here q stands for light quarks and g denotes glu- 
ons. For evaluating r scatt for heavy quarks (Q) the 
pQCD processes are taken from [2g]. The infrared 
divergence appearing in case of massless particle ex- 
change in the i-channel has been shielded by Debye 
mass. 

The expansion time scale can be defined as: 



de(r, r) 



where e(r, r) is the energy density, r and r are proper 
time and the radial co-ordinate respectively. e(r, r) 
is calculated by solving the hydrodynamical equa- 
tion: 



d^ v = 



(4) 



with the assumption of boost invariance along lon- 
gitudinal direction [2{| and cylindrical symmetry of 
the system [3(J. In Eq. Q, = (e + P)u"u" - 
g^"P is the energy momentum tensor, P is the pres- 
sure, m m denotes four velocity and g^ u stands for 
metric tensor. We consider a net baryon free QGP 
here, therefore the baryonic chemical potential (/j-b) 
is zero. 

The expansion rates for RHIC and LHC energies 
have been calculated using the initial conditions, 
T l = 400 MeV, n = 0.2 fm for RHIC which gives 
dN/dy ~ 1100 [3 and T, = 700 MeV, n = 0.08 
fm for LHC giving dN/dy = 2100 The initial 

radial velocity has been taken as zero for both the 
cases. Two sets (SET-I and SET-II) of equation of 
state (EoS) have been used to study the sensitivity 
of the results on EoS. 

SET-I: In a first order phase transition scenario - 
we use the bag model EOS for the QGP phase and 
for the hadronic phase all the resonances with mass 
< 2.5 GeV have been considered H. and SET- 
II: The EOS is taken from lattice QCD calculations 
performed by the MILC collaboration [33|. 
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Expansion (EOS-II) 
free streaming 
Scattering (gluon) 
Scattering (light quarks) 
Scattering (charm) 




e(r, r) dr 



(3) 
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FIG. 1: Expansion time vs scattering time (calculated 
with pQCD process) in RHIC energy. Here the expan- 
sion time scale has been calculated at r — 1 fm. 



In Fig [T] the scattering time scale is contrasted 
with the expansion time scale for two types of EOS 
mentioned above. For the sake of comparison the ex- 
pansion rate for the extreme case of free streaming 
is also displayed. The scattering rates are evaluated 
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FIG. 2: Same as Fig. \T\ with the pQCD cross section 
enhanced by a factor 2 



with pQCD cross sections. The condition for equi- 
libration in Eq. [T] indicates that the gluons remains 
close to equilibrium, however the charm and bottom 
(not shown in the figure) quarks remain out of equi- 
librium during the entire evolution history. 

However, as mentioned in the introduction the 
analysis of the experimental data within the ambit of 
relativistic hydrodynamics suggest that the matter 
formed in Au+Au collisions at RHIC achieve ther- 
malization. One possible reasons for the thermaliza- 
tion to occur is that the partons interact strongly 
after their formation in the heavy ion collisions. It 
is argued in [34| that the onset of thermalization in 
the system formed in heavy ion collisions at rela- 
tivistic energies can not be achieved without non- 
perturbative effects. It has also been shown in [35| 
that a large enhancement of the pQCD cross sec- 
tion is required for the reproduction of experimental 
data on elliptic flow at RHIC energies. Therefore, 
the pQCD cross sections used to derive the results 
shown in Fig. [T] should include non-perturbative ef- 
fects. To implement this we enhance the pQCD cross 
sections by a factor of 2. The resulting scattering 
time is compared with the expansion time in Fig. O 
It is observed that the gluons are kept in equilibrium 
throughout the evolution, light quarks are closer to 
the equilibrium as compared to the heavy flavours. 

In Figs. OH] the results for LHC are displayed for 
the two time scales mentioned above for pQCD and 
enhanced cross sections. The expansion becomes 
faster at LHC than RHIC because of the higher in- 
ternal pressure. As a consequence, it is interesting to 
note that the thermalization scenario at LHC does 
not differ drastically from RHIC. 
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FIG. 3: Same as Fig. Q] for LHC energy. 
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FIG. 4: Same as Fig. [2] for LHC energy. 



III. NON-EQUILIBRIUM PROCESS 



It is argued in [19( that the relaxation time for 
heavy quarks is larger than the corresponding quan- 
tities to light partons by a factor of M/T, where 
M is the mass of the heavy flavour and T is the 
temperature. In the present work we have also seen 
that the heavy flavours do not maintain equilibra- 
tion throughout the evolution scenario, but the glu- 
ons are close to equilibrium. Therefore, we treat this 
problem as an interaction between equilibrium and 
non-equilibrium degrees of freedom and FP equation 
provides an appropriate framework for such studies. 
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The Boltzmann transport equation describing a 
non-equilibrium statistical system reads: 



d p d _ d 
h — h F. — 

dt E dx dp 



f(x,p,t) 



df 
dt 



(5) 



col. 



The assumption of uniformity in the plasma and ab- 
sence of any external force leads to 



df 
dt 



df 



dt 



(6) 



The collision term on the right hand side of the above 
equation can be approximated as (see [Tt], for 
details): 



\df 


d 


.dt_ 


coi d Pi 



_d_ 

dpt ' 



Ai^f + irlBiMf] 



(7) 



where we have defined the kernels 
A, = I d 3 puj(p, k)ki 



B i: j = J d 3 pui(p, k)kikj 



and the function ui(p, k) is given by 
f d 3 q 

u(p, k) = 9j J j^fM^ P , q -^p-k, q +k 



(8) 



(9) 



where fj is the phase space distribution for the par- 
ticle j, Vij is the relative velocity between the two 
collision partners, a denotes the cross section and 
gj is the statistical degeneracy. The co-efficients in 
the first two terms of the expansion in Eq. [7] are 
comparable in magnitude because the averaging of 
ki involves greater cancellation than the averaging 
of the quadratic term kikj. The higher power of ki's 
are smaller (37j . 

With these approximations the Boltzmann equa- 
tion reduces to a non-linear integro-diffcrential equa- 
tion known as Landau kinetic equation: 



dl = _d_ 

dt dpi 



M P )f+^[B tJ (p)f] 



(10) 



The nonlinearity is caused due to the appearance 
of / in Aj and Bij through w(p,k). It arises from 
the simple fact that we are studying a collision pro- 
cess which involves two particles - it should, there- 
fore, depend on the states of the two participating 
particles in the collision process and hence on the 
product of the two. Considerable simplicity may be 
achieved by replacing the distribution functions of 
the collision partners of the test particle by their 
equilibrium Fermi-Dirac or Bose-Einstein distribu- 
tions (depending on the statistical nature) in the 
expressions of A t and . Then Eq. [TU] reduces to a 



linear partial differential equation - usually referred 
to as the Fokker-Planck equation[38[ describing the 
motion of a particle which is out of thermal equi- 
librium with the particles in a thermal bath. The 
quantities Ai and B^ are related to the usual drag 
and diffusion coefficients and we denote them by 
and Dij respectively (i.e. these quantities can be 
obtained from the expressions for Aj and B^ by re- 
placing the distribution functions by their thermal 
counterparts) : 



dl = d_ 
dt dpi 



(11) 



We evaluate the value of the ji and Dij for the reac- 
tion: gQ — > gQ and qQ — > qQ for both zero and non- 
zero quark chemical potential (p, = /is/3). In Fig. [5] 
we depict the variation of the drag coefficients as a 
function of the transverse momentum of the charm 
and bottom quarks at a temperature, T — 200 MeV. 
The momentum dependence is weak. For non-zero 
quark chemical potential the value of the drag in- 
creases, however, the nature of the variations remain 
same. In Fig. [6] the temperature variation of the 
drag co-efficient is plotted for both zero and non-zero 
quark chemical potential. Qualitatively, the inverse 
of the drag co-efficient gives the magnitude of the 
relaxation time. Therefore, the present results indi- 
cate that a system with fixed temperature achieves 
equilibrium faster for non-zero /x. In Fig. [7] the dif- 
fusion coefficients are plotted as a function of pt for 
T = 200 MeV. The diffusion co-efficient for non-zero 
\x is larger as compared to the case of vanishing /i. 
The same quantity is displayed in Fig. [5] as a func- 
tion of temperature. In the present work we confine 
to /i = 0. Recently the heavy quark momentum dif- 
fusion co-efficient has been computed [39] at next to 
leading order within the the ambit of hard thermal 
loop approximations. For T ~ 400 MeV our mo- 
mentum averaged pQCD value of the diffusion co- 
efficient is comparable to the value obtained in [39[ 
in the leading order approximation for the same set 
of inputs (e.g. strong coupling constant, number of 
flavours etc). 

The inverse of the drag co-efficients gives an esti- 
mate of the thermalization time scale. Results ob- 
tained in the present work indicate that the heavy 
quarks are unlikely to attain thermalization at RHIC 
and LHC energies [40| . 

The total amount of energy dissipated by a parton 
depends on the path length it traverses through the 
plasma. Each parton traverse different path length 
which depends on the geometry of the system and on 
the point where its is produced. The probability that 
a parton is created at a point (r, <f>) in the plasma 
depends on the number of binary collisions at that 
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point which can be taken as [211 ]: 

P ( r >*) = ^ (1 -^ )(9(iZ - r) (12) 

where R is the nuclear radius. A parton created at 
(r, 0) in the transverse plane propagate a distance 
i = \J R 2 — r 2 sin 2 (f) — rcos(f> in the medium. In 
the present work we adopt the following averaging 
procedure for the transport coefficients. For the drag 
coefficient (7): 



r = / rdrd(j>P(r, < 



L/v 



(13) 



where v is the velocity of the propagating partons. 
The quantity T appears in the solution of the FP 
equation(see [2(| for details). Similar averaging has 
been done for the expression involving diffusion co- 
efficients to take into account the geometry of the 
system. 

Using the drag and diffusion co-efficients as in- 
puts we solve the to FP equation with the following 
parametrization of the initial momentum distribu- 
tion of the charm quarks generated in p-p collisions 
at a/s = 200 GeV: 



d 2 N c 
dp\ 



= C 



frr + A) 2 



(14) 



where A=0.5 GeV, B=6.6 GeV, a = 21 and C=0.845 
GeV -4 . We do not elaborate here on the procedure 
of solving the FP equation as this has been discussed 
in detail in ref. [2Qj. The corresponding initial dis- 
tributions for bottom quarks can be obtained from 
the results obtained in Ref. [4l| for pp collisions at 
sfs = 200 GeV. Obtaining the solution of the FP 
equation for the heavy (charm and bottom) quarks 
we convolute it with the fragmentation functions of 
the heavy quarks to obtain the px distribution of 
the D and B mesons. The following three sets of 
fragmentation functions have been used to check the 
sensitivity of the results: 
1. SET-I [H 



1 hm 2 

f^^^T^^-zTeM--^) (15) 

where tuq is the mass of the charm quark, tq =1, 
a = 5, b = 1 and m\ = toq + p\. It has been 
explicitly checked that the Raa is not very sensitive 
to the values of a and b. 
2. SET-II [H 

f(z)ocz a (l-z) (16) 
where a = — 1 for charm quark and it is 9 for bottom 
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FIG. 5: Variation of drag coefficient with pr for T = 200 
MeV 
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FIG. 6: Variation of drag coefficient with temperature. 



3. SET-III [4J 

/(*)« 



z z- 



I _ -6£-]21 
z 1-zJ J 



(17) 



for charm quark e c = 0.05 and for bottom quark 
e b = (m c /m b ) 2 e c . 

Recently, the pr spectra of D mesons has been 
measured by the STAR collaboration [45[ in Au+Au 
collisions at ^/snn — 200 GeV. The pr spectrum of 
hadrons can be written as I4a| 



dN 



ft 



d 2 prdy (27r) 3 

—pTCOstfidr) f (m m p m ) 



Trd(pdrj{mTCOsh{ri — y)dr 

(18) 



rj is the space time rapidity, p M is the four momentum 
and w/ 1 = 7(1, 0) is the hydrodynamic four velocity, 
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FIG. 7: Variation of diffusion coefficient with pp for T 
200 MeV. 



FIG. 9: Invariant momentum distribution of the D- 
meson as a function of pp. 
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FIG. 8: Variation of diffusion coefficient with Tempera- 
ture. 



u^p^ is the energy of the hadrons in the co-moving 
frame of the plasma and /{u^p^) is the momentum 
space distribution. In the spirit of the blast wave 
method we can write Eq. [18] as j4?| : 



dN 



fJ 



rrd(j)dri 



d 2 pTdy (2tt) 3 

m T cosh(r] - y)f(u tl p tl )dr 

Taking the surface velocity profile as: 



/ r \ n 

m = A (-) 



(19) 



(20) 



and choosing n=l, the pt spectra of D mesons is 
evaluated. 

Before comparing the data with non-equilibrium 
momentum distribution we analyze the data within 



the ambit of the blast wave method (47J assuming a 
equilibrium distribution for the D-meson. The val- 
ues of the blast wave parameters i.e. the radial flow 
velocity at the surface, (3 S and the freeze-out temper- 
ature, Tp are 0.6 and 0.170 GeV respectively. The 
data is reproduced well (FiglHJ). The value of Tp is 
close to T c , which indicate that the D-mesons (even 
if the charm is in equilibrium in the partonic phase) 
can not maintain it in the hadronic phase. This is 
reasonable because of the low interaction cross sec- 
tions of the D mesons with other hadrons. Next 
we replace the equilibrium distribution in Eq. [TO] by 
the solution of FP equation appropriately boosted 
by the radial velocity. The results are displayed 
in Fig. [S] The data is reproduced well for all the 
three sets of fragmentation functions mentioned be- 
fore. The value of the freeze-out temperature is 170 
MeV and the flow velocity at the surface is 0.45, 
0.35 and 0.4 for SET-I, SET-II and SET-III frag- 
mentation functions respectively. The values of f3 s is 
lower here than the equilibrium case for all the frag- 
mentation function. It is interesting to note that at 
low pr{< 0.5 GeV) domain the results for equilib- 
rium distribution substantially differ from the non- 
equilibrium distribution for all the three sets of frag- 
mentation functions. Therefore, measurements of 
the heavy meson spectra at low pp domain will be 
very useful to distinguish between the equilibrium 
and the non-equilibrium scenarios. The two scenar- 
ios also give different kind of variation at large pp. 
The pp integrated quantity, i.e. the D meson mul- 
tiplicity may also be useful to understand the differ- 
ence between the equilibrium and non-equilibrium 
scenarios. 
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FIG. 10: Nuclear suppression factor, Raa as function of 

PT- 



FIG. 12: Same as Fig. 1111 for fragmentation functions of 
SET-II and III. 
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FIG. 11: Same as Fig. [10] with enhancement of the cross 
section by a factor of 2. 



IV. NON-PHOTONIC SINGLE ELECTRON 
FROM HEAVY FLAVOURS. 

The STAR [H and the PHENIX [49] collabora- 
tions have measured the non-photonic single electron 
inclusive pt spectra recently both for Au+Au and 
p+p collisions at ^/snn — 200 GeV. The quantity 



Raa(j>i 



dN e Au+Au 

d'^pTcly 

N x dN" P+P 
coll d' 2 pTdy 



(21) 



called the nuclear suppression factor, will be unity 
in the absence of any medium. However, the exper- 
imental data from both the collaborations [H, |4^| 
shows substantial suppression (Raa < 1) for px > 2 



GeV indicating the interaction of the plasma parti- 
cles with the charm and bottom quarks from which 
electrons are originated through the process: c(b) 
(hadronization) — ► D(B)(dccay) — ► e + X. The 
loss of energy of high px heavy quarks propagating 
through the medium created in Au+Au collisions 
causes a depletion of high pt electrons. 

To evaluate Raa theoretically, the solution of 
the FP equation for the charm and bottom quarks 
should be convoluted by the fragmentation func- 
tions to obtain the pr distribution of the D and B 
mesons which subsequently decay through the pro- 
cesses: D — > Xev and B —t Xev. Similar formalism 
has been used in [5(| to study the evolution of light 
quark momentum distributions. The resulting elec- 
tron spectra from the decays of D and B mesons can 
be obtained as follows [HJ, [H, [HJ : 



dN e 
prdpT 



where 



dN D 

d<lT — j — F (pt , <Zt ) (22) 
qrdqr 



F(p T , q T )=oj f ,f PT-qx) g (p T .q T /M) (23) 
J 2p T p T .qT 

where M is the mass of the heavy mesons (D or B), 
uj = 96(1 - 8m 2 + 8m 6 - to 8 - UrnHnm 2 )- 1 M~ 6 
(m = mx/M) and g(E e ) is given by 

g{Ee = El(M 2 — M\ — 2ME e ) 2 

related to the rest frame spectrum for the decay 
D(B) — > Xev by the following relation [5l| 



1 dY 

TdEZ = " 9{Ee) - 



(25) 
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We evaluate the electron spectra from the decays 
of heavy mesons originating from the fragmentation 
of the heavy quarks propagating through the QGP 
medium formed in heavy ion collisions. Similarly 
the electron spectrum from the p-p collisions can be 
obtained from the charm and bottom quark distribu- 
tion which goes as initial conditions to the solution 
of FP equation. The ratio of these two quantities 
gives the nuclear suppression Raa- For a static sys- 
tem the temperature dependence of the drag and 
diffusion co-efficients of the heavy quarks enter via 
the thermal distributions of light quarks and gluons 
through which it is propagating. However, in the 
present scenario the variation of temperature with 
time is governed by the equation of state or veloc- 
ity of sound of the thermalized system undergoing 
hydrodynamic expansion. In such a situation the 
quantities like T (Eq. [T3]) and hence Raa becomes 
sensitive to velocity of sound in the medium. 

The results for Raa is displayed in Fig. [TDJ The 
theoretical results are obtained for fragmentation 
function of SET-I 0. The velocity of sound for 
the QGP phase is taken as c s = 1/vi corresponds 
to the equation of state, p = e/4. The results failed 
to describe the data in this case. Next we generate 
Raa by changing the value of c s to l/\/5 and keeping 
all the other parameters fixed. The resulting spectra 
describes the data reasonably well. Lower value of c s 
makes the expansion of the plasma slower enabling 
the propagating heavy quarks to spend more time to 
interact in the medium and hence lose more energy 
before exiting from the plasma which results in less 
particle production at high px- Further lowering of 
c s gives further suppressions. 

However, as mentioned earlier the non- 
perturbative effects are important for the interaction 
of the heavy quarks with the plasma. Therefore, we 
enhance the cross section by a factor of 2 and find 
that the experimental results can also be explained 
by taking an EoS P = e/4 ( Fi EflTj) and keeping 
other quantities like fragmentation functions etc 
unchanged. The ideal gas EoS P = e/3 can not 
reproduce the data even if the cross section is 
enhanced by a factor of 2. With c s = 1/4 and 
enhanced cross section (by a factor 2) the data can 
also be described with for fragmentation functions 
of set-II and set-Ill Fig. [T2] Several mechanisms 
like inclusions of non-perturbative contributions 
from the quasi-hadronic bound state [54| . 3-body 
scattering effects [55| and employment of run ning 
coupling constants and realistic Debye mass [57| 
have been proposed to improve the description of 
the experimental data. It is demonstrated here that 
the EoS of the medium and the non-perturbative 
effects play a crucial role in determining the nuclear 
suppression factor. 



V. SUMMARY AND CONCLUSIONS 

The transverse momentum spectra of the D and 
B mesons have been studied within the ambit of 
Fokker-Planck equation where the charm and bot- 
tom quarks are executing Brownian motion in the 
heat bath of light quarks and gluons. We have eval- 
uated the drag and diffusion co-efficients both for 
zero and non-zero quark chemical potential. Re- 
sults for non-zero baryonic chemical potential will be 
very useful for studying physics at low energy RHIC 
run (56j . The results are compared with experi- 
mental data measured by STAR collaboration. It 
is found that the present experimental data can not 
distinguish between the pt spectra obtained from 
the equilibrium and non-equilibrium charm distri- 
butions. Data at lower px may play a crucial role 
in making the distinction between the two. Since 
the results for equilibrium and non-equilibrium sce- 
narios differ, the px integrated quantity i.e. the D 
meson multiplicity may also be very useful to put 
constraints on the model. The nuclear suppression 
factor for the measured non-photonic single electron 
spectra resulting from the semileptonic decays of 
hadrons containing heavy flavours have been eval- 
uated using the present formalism. The experimen- 
tal data on nuclear suppression factor of the non- 
photonic electrons can be reproduced within this for- 
malism if the expansion of the bulk matter is gov- 
erned by a equation of state p = e/4 and the par- 
tonic cross sections are taken as 2 x (j p qcd- Three 
different kinds of fragmentation functions for the 
charm and bottom quarks hadronizing to D and B 
mesons respectively have been used and found that 
c s ~ 1/a/4 can describe the data reasonably well. 
The data can not be reproduced with c s ~ l/\/3 
even after enhancing the cross section by a factor 
of 2. The loss of energy by the heavy quarks due 
to radiative process may be suppressed due to dead 
cone effects. In the present work the radiative loss 
is neglected. The FP equation needs to be modi- 
fied to include the radiative loss [El, [59|, [6(| [U, HI1 
(see [63j j for a review), work in this direction is in 
progress (56j . 

The calculations may be improved by making the 
space time evolution picture more rigorous as follow. 
In the absence of any external force the evolution of 
the heavy quark phase space distribution is governed 
by the equation: 

(| + vp- V r ) /(p, r, t) = C[f(p, r, t)] (26) 

As mentioned before the FP equation can be ob- 
tained from the above equation by linearizing the 
collision term, C[f(p,r,t)]. To take into account the 
energy loss of the heavy quarks in the thermal bath 
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the ideal hydrodynamic equation needs to be modi- 
fied as: 

d^T^ = J v (27) 

containing a source term J u corresponding to the 
energy-momentum deposited in the thermal system 
along the trajectory of the heavy quark, which may 
be taken as J v ~ dp u /dr [IH (jp u is the four mo- 
mentum vector). Eq. (f27f should solved for T(r,t) 
with appropriate equation of state which can be used 
to obtain the surface of hadronization by setting 



T(r c ,t c ) = T c . Subsequently the solution of Eg. |2T)I 
/(r c ,i c ) f° r the heavy quark on this surface should 
convoluted with the fragmentation function to ob- 
tain the B and D distribution. 
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